看了一下vasp加入范式力修正的东西,计算的结果感觉有点错误,不知道错在哪里请各位帮助。
本人计算金属表面吸附气体小分子
修正前的吸附能为-0.02eV
修正后的吸附能为-0.3eV
吸附能的公式:
吸附能E=E(吸附体系能量)-E(分子体系能量)-E(没有吸附的表面能)
如果E(没有吸附的表面能)取修正前的话,E(吸附体系能量用修正后的)就会得到
-4.4eV的吸附能,比化学吸附的分子还强。
如果E(没有吸附的表面能)取修正前的话,吸附能就是-0.3eV。
问题是:范式力的修正不是针对分子间的长程作用么?为什么金属表面加入范式修正后能量会产生变化呢?在算吸附能的后面减去的表面能量是修正前的还是修正后的。如果范式力的作用使表面能量差这么多,那以后的化学吸附是不是也要考虑啊?
下面为我取的范式力的参数:
LVDW=.Ture.
VDW_RADIUS=10 (因为真空是12A,为了避免考虑进去两层间的作用所以将考虑范围缩小到了10A. 默认是30A都测试过能量差0.1eV)
VDW_SCALING=0.75
VDW_D=20.0
VDWC6=1.75 0.14 24.67 10.80
VDW_R0=1.452 1.001 1.639 1.562
望高手赐教! 返回小木虫查看更多
|